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ABSTRACT 

O 

' We discuss non-thermal emission mechanism of the Crab-hke pulsars with 

^ both a two-dimensional electrodynamical study and a three-dimensional model. 

■ We investigate the emission process in the outer gap accelerator. In the two- 
dimensional electrodynamical study, we solve the Poisson equation of the accel- 
erating electric field in the outer gap and the equation of motion of the primary 

■ particles with the synchrotron and the curvature radiation process and the pair- 
j_| ' creation process. We show a solved gap structure which produces a consistent 
^ gamma-ray spectrum with EGRET observation. Based on the two-dimensional 

model, we conduct a three-dimensional emission model to calculate the syn- 
chrotron and the inverse-Compton processes of the secondary pairs produced 
outside the outer gap. We calculate the pulse profiles, the phase-resolved spec- 
tra and the polarization characteristics in optical to 7-ray bands to compare the 
observation of the Crab pulsar and PSR B0540-69. For the Crab pulsar, we find 
that the outer gap geometry extending from near the stellar surface to near the 
light cylinder produces a complex morphology change of the pulse profiles as a 
function of the photon energy. This predicted morphology change is quite similar 
with that of the observations. The calculated phase-resolved spectra are consis- 
tent with the data through optical to the 7-ray bands. We demonstrate that 
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the 10~20 % of the polarization degree in the optical emissions from the Crab 
pulsar and the Vela pulsar are explained by the synchrotron emissions with the 
particle gyration motion. For PSR B0540-69, the observed pulse profile with a 
single broad pulse is reproduced with a thicker emission region and a smaller in- 
clination angle between the rotational axis and the magnetic axis than the Crab 
pulsar. 

Subject headings: optical-X ray-gamma rays:theory-pulsars:Crab like -radiation 
mechanisms:non-thermal 



1. Introduction 

The observed strong 7-ray emissions from the seven young pulsars (Thompson 2003) 
show that electrons and positrons are accelerated up to ultra-relativistic regime in the pulsar 
magnetosphere. The Crab pulsar (PSR B0531-I-21), which is one of the brightest and the 
youngest 7-ray emitting pulsar, shows the non-thermal emission properties in optical to 7- 
ray bands. The observed spectrum of the pulsed photons emitted from the Crab pulsar 
extends continuously from optical to 7-ray bands with the spectral index a^, defined as 
Ii, oc z/^"", varying from ~ in optical wavelengths, a^, ~ 0.5 in X-ray bands, to cij, ~ 2 
in 7-ray bands. The pulse profile has two peaks in a single period, and the positions of 
the pulse peaks across the wide energy range arc approximately all in phase (Kuiper et 
al. 2001). Interestingly, the pulse profile morphology changes significantly as a function of 
the photon energy. The first peak (denoted Peak 1 in the following) dominates in optical 
wavelengths. However, the second peak (Peak 2) becomes more and more pronounced for 
increasing energies and eventually the Peak 2 emission dominates in soft 7-ray bands. Above 
10 MeV photon energy. Peak 1 again dominates Peak 2. The electromagnetic spectrum of 
the non-thermal emissions also changes with pulse phases. In the future, the phase-resolved 
spectra above 10 McV will be measured with a sensitivity better than that of the Energetic 
Gamma-Ray Experiment Telescope on board the Compton Gamma-ray Observatory by, for 
example, GLAST LAT. These observed detail properties for the pulse profiles and the phase- 
resolved spectra will be useful to discriminate the proposed emission models. 

In addition to the pulse profiles, Kanbach et al.(2005) measured the polarization char- 
acteristics of the pulsed photons from the Crab pulsar in the optical wavelengths. The 
observation revealed that the degree of the polarization at each pulse peak is lower than 
10 % and a large swing of the position angle of the electric- vector of the radiation appears at 
each pulse peak. The polarization measurements provide two additional observed properties, 
namely, the degree and the position angle of the polarization. In the future, the polarization 
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of X-ray and soft 7-ray emissions from the pulsar will probably be able to be measured 
by ongoing projects such as PoGO (Kataoka et al. 2005) and NCT (Chang et al. 2007) 
projects. Therefore, a theoretical study, particularly on the polarization characteristics, is 
not only desired, but also timely for the understanding the non-thermal emission process in 
the pulsar magnetospheres. 

The polar cap accelerator (Ruderman and Sutherland 1975; Daugherty and Harding 
1996) and the outer magnetospheric accelerator, the so called outer gap model (Cheng et 
al. 1986a,b; Romani 1996), were proposed as the possible acceleration sites in the pulsar 
magnetospheres. The traditional polar cap model assumes an acceleration region expanding 
several stellar radii from the stellar surface around the magnetic pole. On the other hand, 
the traditional outer gap assumes an acceleration region extending beyond the null surface 
of the Goldreich-Juhan charge density at the outer magnetosphere. The Goldreich- Julian 
charge density is given by poj ~ — • B/27rc (Goldreich and Julian 1969) with Q, being the 
rotational frequency of the star, B the magnetic field, and c the speed of light. Both models 
assume the particle acceleration by an electric field parallel to the magnetic field line. In 
the pulsar magnetosphere, the accelerating electric field arises in the region where the local 
charge density differs from the Goldreich- Julian charge density. 

The slot gap model (Muslimov & Harding 2004), which is an extended polar cap model, 
predicts that the acceleration region extends up to near the light cylinder around the last- 
open field lines because the pair-formation front, which screens the accelerating electric field, 
occurs at higher altitude around there. Two-dimensional electrodynamical studies (Takata 
et al. 2004, 2006; Hirotani 2006) suggested that the inner boundary of the outer gap locates 
near (or at) the stellar surface because of the current through the outer gap. Although the 
recent polar-slot gap and outer gap models both predict similar geometry of the acceleration 
region, an important difference between the two models is the electric field configuration in 
the accelerator. For the slot gap accelerator, the electric field is stronger nearer the stellar 
surface and smaller at higher altitude. On the other hand, the outer gap model predicts a 
stronger electric field beyond the null surface and a smaller one below the null surface due 
to the screening effect of electron and positron pairs. This difference in the electric field 
configuration, and the resultant difference in the acceleration and the emission structures 
will appear as a difference in the predicted pulse profiles, the phase-resolved spectra and the 
polarization characteristics, which can be examined by a three-dimensional model. 

Within the framework of the traditional outer gap model, Romani & Yadigaroglu (1995) 
considered a three-dimensional geometry and explained the general features of the observed 
pulse profile such as two-peaks in a single period. Subsequently, Cheng et al. (2000, here- 
after CRZOO) developed the three-dimensional outer gap model, in which the gap is sustained 
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self-consistently by the pair-creation process between the primary photons emitted via the 
curvature process and the surface X-ray photons. CRZOO calculated the phase-resolved 
spectrum in 7-ray region for the Crab pulsar with the synchrotron radiation and the inverse 
Compton scattering of the electron and positron pairs produced outside the gap. Zhang & 
Cheng (2002) reconsidered the CRZOO model to examine the phase-resolved spectra in X-ray 
region. However, it has been difficult to explain the detail features of the observed pulse 
profiles and phase-resolved spectra with the traditional model. Furthermore, the previous 
studies have not discussed the complex features of the observed pulse profiles from optical 
to 7-ray bands. Recently, Takata et al. (2007) explained the observed polarization charac- 
teristics in the optical wavelengths (Kanbach et al 2005) with the new outer gap geometry. 
Jia et al. (2007) examined the phase-resolved spectra by taking account of the emissions 
below null charge surface. However, these studies also did not consider the pulse profile, 
the phase-resolved spectra and the polarization characteristics in optical to 7-ray bands, 
simultaneously. 

In this paper, we study the emission process of the Crab-like pulsars with the outer 
gap accelerator model from both a two-dimensional electrodynamical model and a three- 
dimensional emission model point of views. In first part (section [2]) of this paper, we will 
summarize the results of the two-dimensional electrodynamical study, in which the outer 
gap structure for the Crab pulsar is solved with the Poisson equation, the particle motion, 
the radiation process and the pair-creation process in meridional plane, following Takata 
et al (2004, 2006) and Hirotani (2006). We will show a result which has a consistent GeV 
spectrum with the observed phase-averaged spectrum of the Crab pulsar. In the second 
part (sections [3] and Hj), we will conduct a three-dimensional outer gap model based on the 
results of the two-dimensional electrodynamical study. In the three-dimensional study, the 
main purpose is to discuss the emission process of optical to 7-ray photons by examining the 
morphology change of the pulse profile as a function of the photon energy and the phase- 
resolved spectra for the Crab pulsar with the outer gap accelerator model. We will predict 
the polarization characteristics through optical to 7-ray bands. We also apply the model 
to a Crab-like pulsar, PSR B0540-69. The Crab pulsar and PSR B0540-69 are sometimes 
called twin pulsars, because their pulsar parameters are very similar to each other. However, 
the observed shapes of pulse profiles are very different to each other. This pair will give an 
unique opportunity to examine the model capability. 

Important differences between present and previous three-dimensional studies are as 
follows. First, we take into account the emissions both below and beyond the null surface as 
the electrodynamical study has predicted, while only the emissions beyond the null surface 
were taken into account in CRZOO. Secondary we discuss the morphology change of the 
pules profile by calculating local emissivity as a function of the photon energy, while the 
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previous studies did not discuss the morphology change because they assumed a constant 
emissivity when the pulse profiles were calculated. We deal the gyration motion of the pairs 
because the gyration motion causes the depolarization for the synchrotron radiation. We 
adopt the rotating dipole field in the observer frame, while the previous studies adopted it 
in the co-rotating frame. Though these effects were considered in Takata et al. (2007), they 
calculated only the synchrotron emission process and presented the phase-averaged spectrum 
below MeV energy. In this paper, we extend the model spectrum up to 7-ray bands by 
computing also the inverse Compton scattering. Finally, we calculate the collision angle of 
the inverse Compton scattering between the pairs and the background synchrotron photons 
by tracing the three-dimensional trajectory of the synchrotron photons, while the isotropic 
distribution of the back ground photons was assumed in the previous studies (CRZOO). The 
collision angle greatly affects to the emissivity of and the polarization characteristics of the 
inverse Compton scattering. By including all these effects, we examine the pulse profiles, 
the phase-resolved spectra and the polarization characteristics in optical to 7-ray bands, 
simultaneously. 



2. Results of Two-dimensional Electrodynamical Model 

In this section, we summarize the results of the two-dimensional electrodynamical study 
for the Crab pulsar. Following Takata et al (2004, 2006) and Hirotani (2006) we calculate 
the spectrum of the synchrotron and curvature radiation processes of the primary particles 
with the electric structure by solving the Poisson equation [V^$ = — 47r(p — pcj)], the 
equation of motion for the particles, the pair-creation process and the radiation process. As 
discussed in Takata el al. (2004, 2006), the electric structure depends on the current and 
the gap size, which are model parameters in their studies. In this section, we show a result, 
which produces a consistent GeV spectrum with the observations. We ignore the effect of 
the gravity which is not important for the dynamics of the outer gap accelerator. We adopt 
static dipole field, while in the later section of the three-dimensional study, we apply the 
rotating dipole field. The obtained electric structure with the static and the rotating dipole 
field did not change very much, because the radial distances to the null charge points, that 
is, to the gap position are similar to each other. For the pair-creation process in the gap, we 
consider the thermal soft-photons coming from the stellar surface. We adopt kT =170 eV 
for the Crab pulsar (Yakovlev & Pethick 2004). The inclination angle is assumed as a = 50°. 

Thick solid line in Figured] shows the solved accelerating electric field along the field line 
locating at 50 % of the trans- field thickness from the lower boundary (last-open field line). 
Here, we assume O.lRic of the gap thickness at the hght cylinder and the outer boundary is 
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putted at near the light cyhnder. We also assume that 5% of the Goldreich- Julian current, 
0.05flB/2n, is injected at the outer boundary. 

The position of the inner boundary is solved with the current, for which about 22 % of the 
Goldreich- Julian current runs through the outer gap in the present case. From Figure [H we 
can see that the inner boundary (r ~ O.lSRic) is inside of the null charge point (r ~ 0.29Ric)- 
As suggested by Takata et al. (2004), the inner boundary is located at the position, on which 
jg+j2—ji ~ Bz/B is satisfies, where jg is non-dimensional current created in the gap, ji and 
j2 are non-dimensional current injected at the inner and the outer boundaries, respectively. 
For example, for no injection currents, ji = j2 = 0, the inner boundary is located at the 
null charge surface, where 5^ = 0, if no current is created inside of the gap (jg = 0) as 
the vacuum case. On the other hands, if jg ~ cos a, where a is the inclination angle, is 
created, the inner boundary is located at the stellar surface on which B^/B ~ cos a is 
satisfied around the magnetic pole. In the present case, the inner boundary is located at the 
position of about 65 % of the radial distance to the null point with the current components 
(jg, Ji, J2) = (0.17, 0, 0.05). 

Figure [2] shows the calculated synchrotron-curvature spectrum and compares with the 
observed phase-averaged spectrum. Sold-line shows spectrum of the intrinsic radiation from 
the outer gap, while the dashed-line represents the appearance spectrum after attenuation 
of the photons via pair-creation process outside of the gap with the soft photon-field emitted 
by the synchrotron process of the secondary pairs. We calculate the initial pitch angle of 
the pairs from the propagating direction of the curvature photons and the magnetic field 
direction at the pair-creation position. For obtaining the luminosity, we assume the gap 
opening angle ~ 250 degree in the azimuthal direction (see section ISTTj) . 

As dashed-line shows the large amount of the curvature photons above 500 MeV are 
converted into the secondary pairs outside of the gap via the pair creation process with 
the X-ray photons from the secondary pairs. We find from Figure [2] that the shape of 
the spectral energy distribution after absorption becomes relatively flat and explains the 
observation above 100 MeV. 

The inverse- Compton process of the primary particles in the gap is a possible mech- 
anism for TeV emissions. However, the present model predicts the TeV flux for the Crab 
pulsar is too low to detect the present Cherenkov telescopes. The soft-photons emitted by 
secondary pairs above the gap may not be able to illuminate the gap due to the curvature 
of the field lines, and only thermal photons from the stellar surface may be scattered by 
the primary particles. In such a case, we found that the intrinsic flux on the Earth be- 
comes ~ 10~^^ erg/cm^s , which is much small compared with the sensitivity of the present 
Cherenkov telescope. On the other hand, some soft-photons emitted by secondary pairs may 
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illuminate the outer gap because of the effects of the pith angle. In such a case, we obtained 
intrinsic TeV flux which is easily detected by present instruments. However, the optical 
depth of the pair-creation for TeV photons is much larger than unity in the magnetosphere, 
and the residual TeV photons are a very few, which is difficult to detect with the present 
instruments. 

From Figure [21 we can see that the synchrotron and curvature radiations of the primary 
particles in the outer-gap does not explain the observed flux below 100 MeV. We consider that 
the secondary pairs created outside gap produce below 100 MeV photons via the synchrotron 
and the inverse-Compton process. Furthermore, the present two-dimensional model can 
compare with only the phase-averaged spectrum. More detailed observation such like the 
pulse profile, the phase-resolved spectra and the polarization require a three-dimensional 
model. Following sections, therefore, we calculate the emission process of the secondary 
pairs and conduct a three-dimensional model. 



3. A Three-Dimensional Emission Model 



In the following, we conduct a three-dimensional emission model. We anticipate that 
the emission direction is coincide with the particle motion in the observer frame. In the 
present paper, we adopt the rotating dipole field in the observer frame while it was assumed 
in the co-rotating frame in the previous studies (Romani & Yadigaroglu 1995; Cheng et al 
2000; Dyks et al 2004). As a results, the magnetic field configurations and resultant the 
morphology of emission pattern in the observer frame are different between the present and 
the previous studies, but the difference becomes to be important only near the light cylinder. 
In the present study, furthermore, we discuss the model in the observer frame only, and we 
do not introduce the co-rotating frame. 



3.1. Electric field 

We have to describe the accelerating electric field into the three-dimensional from. Based 
on the result (Figure [T]) of the two-dimensional electrodynamical study, we adopt the fol- 
lowing three-dimensional form. First, we use the vacuum solution obtained by Cheng et al. 
(1986a), 

^ OBMl, ,1) 

cs[r) 

beyond the null charge surface, where s{r) is the curvature radius of the magnetic field line 
and /(r) is the fractional gap thickness. We can calculate the electric field at each point 
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having a three-dimensional radial distance, r. Below null surface, we assume 50 % of the 
strength of the electric field at the null point given by equation ([1]). The dashed-line in 
Figure [T] shows the electric field strength of the approximation form in the meridional plane. 
We can see that the typical strength of the electric field by the two-dimensional electro- 
dynamical study (solid-line) in the meridional plane is in general described by the present 
simple form (dashed-line). In fact, as long as the gap is geometrically thin in the trans-field 
direction and the magnitude of the current is smaller than the Goldreich- Julian value, the 
vacuum solution beyond the null charge surface approximately describes the typical strength 
of the accelerating field in the meridional plane. Now, we assume that this simple form can 
describe also the typical strength of the three-dimensional distribution of the accelerating 
electric field. 

At each point, the maximum Lorentz factor of the particles are determined by the 
force balance between the acceleration by the electric and curvature radiation back re- 
action, Tp{r) = [3s2(r)E||/2e]l/^ where Q2 = fi/lOOs"^ . The primary particle emits 
high-energy photons as the curvature radiation process, whose typical energy is Ecurvij) = 
3hrp{r)c/4:ns{r), and the local power of the curvature radiation is given by Icurv = eE\\c. 

It is important to estimate the polar cap opening angle of the active region of the outer 
gap accelerator. When we consider an open-field line through the outer gap, the pair-creation 
process between the primary curvature photons and the surface X-rays mainly occurs near 
and below the null charge surface. Therefore, we may be able to relate the opening angle with 
the pair-creation mean free path, which is estimated as Z(r) ~ [2s{r)f{r)Ric\ ~ 2f^^'^(Ric/2)r, 
at the null surface. In the present paper, we constrain the width of the polar cap angle of the 
active gap by the condition that the mean-free path at the null charge point on the magnetic 
field line becomes shorter than the light radius. This condition produces the width of the 
polar cap angle of ~ 250°. 



3.2. Distribution and motion of the secondary pairs 

As we demonstrated in section [21 a significant amount of the curvature photons above 
~500 MeV convert into secondary pairs outside the gap via the photon-photon pair-creation 
process with the soft-photons emitted by the synchrotron radiation of the secondary pairs. 
From Figure [21 furthermore, we can read that the spectrum of the intrinsic emissions has the 
photon index of about —1. Therefore, we may approximately describe the local curvature 
spectrum with 

Fcurv ~ ^curvj^GJ / EcurvE^, wherc j represent the current in units of the 
Goldreich- Julian value. Using the steady loss equation, d[Eedn/dEe]/dEe = Q{Ee), we 
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obtain the distribution of the secondary pairs as 

drie ^ ( lcurvjnGjMEcurv/'2Ee)/EeEcurv for 500 MeV < 2Ee < E^urv 

dE^ ^ \ h^rvjnGjHEcurv/^QQ MeV)/E,E,^rv for 2E, < 500 MeV ^ ^ 

where i?e = 2e^i?^(r) siia^ 9p{r)Tl/3mlc^ is the energy loss rate of the synchrotron radiation 
of the secondary pairs, 6p is the pitch angle, and Fg is the Lorentz factor of the pairs. 
The local pitch angle will be expressed as sin 6'p(r) oc \j2f{r)RL/s{r). Because /(r) = 
f{Ri^{r/RicY^'^ and s(r) ~ VrRj2 are satisfied for the dipole field, we may relate with 
sin 6'p(r) = (r / RicY^'^ sin 9p{Ric) between the pitch angle of the local point and the light 
cylinder along the field lines. Outside the gap, the pairs loose most of their energy via the 
synchrotron process. Because the synchrotron loss rate E^ is proportional to square of the 
particle energy, the power law index of the distribution given by equation ([2]) becomes p ~ 2, 
which produces a synchrotron spectrum with the spectral index of ~ 0.5. 

In the observer frame, we may describe the particle motion outside the gap with 

P = Po cos 9pb + Po sin 9pb^ + Pco^^, (3) 

where the first term and the second term in the right hand side represent, respectively, 
the particle motion parallel to the magnetic field and gyration motion, and the third term 
represents the co-rotational motion with the non-dimensional velocity f3co = ^/Ric, where 
w is the axial distance. The vector b is the unit vector along the field line and 6^ represents 
the unit vector perpendicular to the magnetic field line, b^ = ±{cosS(j)K + sin6(j)K x b), 
where the sign + (or — ) corresponds to gyration of the positrons (or electrons), K = 
(6 ■ V)b/|(6 ■ V)6| is the unit vector of the curvature of the magnetic field line, and 50 
represents the phase of gyration around the magnetic field. Because the pairs have an 
ultra-relativistic speed, we determine the value of the coefficient /3o from the condition that 
1/3 1 = 1. We anticipate that the photons are emitted in the direction of the particle motion 
of equation ([3]). 

We note that the synchrotron radiation after collecting of the photons is greatly depo- 
larized due to the gyration motion of the pairs, although the intrinsic radiation is highly 
polarized. Therefore, the observed small polarization degree ~ 10 % at the optical bands 
for not only the Crab pulsar, but also for the Vela pulsar (Mignami et al. 2007) are easily 
reproduced by the synchrotron emission model (Takata et al. 2007). 

We assume that the emission region of the secondary pairs extends just above the outer 
gap with thickness of the mean free path of the pair-creation A ~ 10'' cm~ O.lRic for the 
Crab pulsar. Some secondary high-energy photons via the inverse Compton scattering may 
convert into the tertiary pairs. The tertiary pairs will be produced above the emission region 
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of the secondary pairs, and its initial Lorentz factor will be smaller than that of the secondary 
pairs. We also take into account the effects of the emissions from the tertiary pairs. 



3.3. Emission process of the secondary pairs and polarization 

We consider that the synchrotron radiation and the inverse Compton scattering of 
the secondary pairs are major emission mechanisms for the observed non-thermal radia- 
tion through optical to 7-ray bands for the Crab-like pulsars. If we estimate the radiation 
powers of the synchrotron radiation and the inverse Compton scattering, we obtain 

Psyn _ , 1 n ^pf^ ^ V ^ y /^sin 



Pic V5- lO^erg/cm^J VlO^GaussJ V 0.1 ) 

with Uph being the energy density of the synchrotron photons. The estimated value will 
explain the observed flux ratio of 1 MeV and 100 MeV emissions of the Crab pulsar. We 
calculate only the outward emissions, because the inward emissions are expected to be much 
fainter than the outward emissions. 

In the calculation, we firstly compute the volume emissivity of the synchrotron radia- 
tion and its emitting direction for each radiating point (section 13.3.11) . Then, we trace the 
propagation of the synchrotron beam to simulate the scattering process by the pairs (sec- 
tion l3.3.2|) . On each scattering point, we calculate the volume emissivity and the polarization 
of the inverse Compton scattering for a specific viewing angle of the observer. We perform 
this procedure for all calculation points to obtain the total radiation for the specific observer. 
This procedure is equivalent with computing the radiation transfer, 

1 = Js(fci,ei) +Ji(fci,ei), (5) 

as 

where I{ki, ei) is the total intensity of the beam propagating in the direction of fci, ei is the 
energy of photons in units of the electron rest mass energy, js{ki, ei) is the volume emissivity 
of the synchrotron radiation, and jj(fci,ei) represents the amount of the scattered photons 
into the direction ki and to the energy ei. We neglect the effects of the absorption, because 
the synchrotron self-absorption is not important above optical photon energy, where we are 
now interested in. Also, we ignore the effects of the scattering off from the direction ki of 
the synchrotron photons, because the scattered photons are tiny amounts of total number 
of the synchrotron photons. 
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3.3.1. Synchrotron radiation 

Assuming that all the synchrotron photons are radiated toward the particle motion 
direction, ki = f3, the volume emissivity of the synchrotron radiation is calculated from 

where Sfl is the solid angle of the radiation, and Fsyn is the photon spectrum described by 

S^/VSM sin^Jr) 



/ 


dne{r) 




dEe 



F{ei/e,yn)dEe, (7) 



where esyn{r) = 3heTl{r)B{r) sin 9p{r)/4:'7imlc^ is the typical photon energy of the pairs in 
units of the electron rest mass energy, Fg represents the Lorentz factor of the secondary pairs 
and F{x) = x K^/3{y)dy with K^/^ being the modified Bessel function of order 5/3. 

When we calculate the polarization of the synchrotron radiation, we anticipate that 
direction of the electric vector of the electro-magnetic wave propagating toward the observer 
is parallel to the projected direction of the acceleration of the particle on the sky, Eem oc 
a—{ki-a)ki, (Blaskiewicz et al. 1991), where the acceleration vector a derived from equation 
([3]) is approximately written by a ~ ±i3oljb sin 9p{~ sin 6(f)K + cos6(t>K x 6), where ub is 
the gyration frequency. We assume that the radiation at each point is linearly polarized 
with degree of Hsyn = + + 7/2), where p is the power law index of the particle 

distribution. Because the observed radiation is consist of the radiations from the different 
particles with the different pitch angle, we assume that the circular polarization will cancel 
out and become zero in the observed radiation. The Stokes parameters Qsyn and Usyn are, 
respectively, calculated from dQsyn{ki,ei)/ds = js{ki,ei) cos2r]s{r) and dUsyn{ki,ei)/ds = 
js{ki, ei) sin 2r]s{r), where //^(r) is the position angle defined by the angle between the electric 
vector of the wave and the projected direction of the rotation axis on the sky, flp = Q — 

(fci ■ rj)fci. 



3.3.2. Inverse Compton scattering 

To simulate the scattering process, we trace the three-dimensional trajectory of the syn- 
chrotron photons. When we trace the trajectory of the synchrotron photons, we define the 
Cartesian coordinate such that z-axis is along the rotation axis and the x-axis is in the merid- 
ional plane. By ignoring bending of the trajectory due to the gravity, the position of the pho- 
tons after traveling distance 5s is x{5s) = (xq + k^^Ss) cos{6s/Ric) + {yo + ky^^Ss) sin{6s/Ric), 
y{Ss) = —{xo + kxg6s) sin{6 s/Ric) + {yo + ky^6s) cos{6s/Ric) and z{6s) = z^ + k^^^^s, where the 
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coordinates (xq, yo, ^o) are the radiating point of the synchrotron photon, and (fc^^y, ky^, k^^) 
represents the emission direction at {xq, yo, Zq). The emission direction of the background 
synchrotron radiation is calculated from equation [31 which takes account the aberration due 
to the corotating motion. 

Because the mean free path of a synchrotron photon of the scattering is much longer than 
the light radius, one can consider that the scattering rate is constant along the path of the 
synchrotron photons in the magnetosphere as the first order approximation. We determine 
the scattering points at regular interval, which is much shorter than the gap size, along the 
path of the synchrotron photons. In the calculation, we first compute the Stokes parameter 
of the Compton process in the electron rest frame, and then we transform it to the observer 
frame. In the following, the prime and 'non'-prime quantities represent the quantities in the 
electron rest frame and the observer frame, respectively. A detail derivation of equations of 
the Stokes parameters ([8]) and (fTOj) are seen in Appendix Rl 

We denote the specific intensity of the synchrotron radiation propagating to the direction 
ko in the observer frame with Jo(fco,eo), where eo represents the energy of the background 
photons in units of the electron rest mass energy. In the electron rest frame, the background 
radiation becomes /Q(A;Q,eQ) = DflQ{ko,eo), where Di = e'^/eo = r~^(l + /5cos6'g) is a 
Doppler factor, f3 is the velocity of the scattering particles in units of the speed of light, 
and Fe = 1/ \/l — (3'^. The polar angle 6q is defined by the angle between the directions 
of the particle motion and of the propagation of the background radiation, which becomes 
^0 ~ 0.1 — 0.3 radian in numerically. For the particles with the Lorentz factor 10^ ~ 10^, 
optical to X-ray photons are mainly scattered. In this photon energy bands, the synchrotron 
photons are distributed with a spectral index of a^, ~ 0.5 because the cut-off energy of the 
synchrotron spectrum is ~ 1 MeV and because the particles are distributed with the index 
p ~ 2 (section [3.21) . Because the synchrotron beam from each position is strongly collimated, 
we approximately describe the background beam, in which the center of the beam is direct to 
the polar angle measured from the electron motion direction and the azimuthal direction 
00, as lo^ko, 6) = Coe"°-^5(6' — 9o)6{(f) — 0o), where Cq is evaluated from equation ([6]). 

We are interested in the inverse Compton scattering with the background synchrotron 
radiation, which is partially polarized with 11^^^ ~ 70%. With unpolarized components of 
the background radiation propagating to the direction ko, the volume emissivity ju and the 
Stoke parameter Qu and Uu of the scattered radiation propagating to the direction ki are 
calculated from (see appendix [XI) 



ju{ki,ei) = dlu{ki,ei)/ds 
dQu{ki,ei)/ds 
dUu{ki, ei)/ds 
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and 



-T + -r — sm w' 

?«cos2C-< sin 2C 
sin 2C + cos 2C 



q'^ = sin^ cos 21]', 
u'^ = sin^ sin 2r]', 



(9) 



where the Stokes parameters are measured from the rotation axis of the pulsar projected 
on the sky, and ( is defined by the angle between the directions of the rotation axis and 
the particle motion projected on the sky (Figures [TTj) . The polar angle 6i represents the 
propagating direction of the scattered photons measured from the particle motion direction, 
Ws is the scattering angle defined by cosw^ = ko ■ ki and the azimuthal angle rj is the angle 
between the orthogonal direction to the scattering plane and the direction of the particle 
motion projected on the sky. 

For the polarized component of the background radiation, the volume emissivity and 
the Stokes parameters are calculated from 



jp{ki,ei) = dlp{ki,ei)/ds 
dQp{ki,ei)/ds } = Usyn 

dUp{ki, ei)/ds 



16^' 



/ dV, 



X 



r4+<^(l -/3cos^i)2(l + /5cos^^)«+i V^c 



drie 
dVp 



^ + ^ - sin^ w' cos^ A' 

q'p cos 2C — u,p sin 2( 
q'pSin 2C + u'p cos2C, 
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and 



I sm^ wl 



[sh 
[sh 



I sm^ w[ 



[1 + cos^ w'^) cos2Ap] cos 21]' — 2 cos sin 2Ap sin 21]', 
[1 + cos^ w'J cos 2Ap] sin 21]' — 2 cos w'^ sin 2Ap cos 2eri', 



where Xp is the angle between the polarization plane of the background radiation and the 
plane of the scattering. By exploring the additive property of the Stokes parameters, the 
total volume emissivity and Stokes parameters are, respectively, given by jj = + jp, 

dQi/ds = dQu/ds + dQp/ds and dUi/ds = dUu/ds + dUp/ds. 

After collecting all photons from the possible points for each rotation phase $ and a view- 
ing angle ^, the degree of the radiation and the position angle of the electric vector of the radi- 
ation are, respectively, calculated from P(^,$,ei) = ^J Q'^{^, ei) + U'^{^, $, ei)//(^, $, ei) 



-14- 



and ^, ei) = 0.5atan[f/(e, ei)/Q(e, d)], where /(e, ei), Q(e, ei) and f/(e, $, ei) 
are the Stokes parameters after collecting photons emitted via both synchrotron radiation 
and the inverse Compton scattering. The position angle xi^^ ^i) is measured anticlockwise 
from the axis of the rotation projected on the sky (Figure [TT]) . 

3.4. Model parameter 

The inclination angle of the pulsars has been constrained by the polarization measure- 
ments of the radio pulsed emissions. However, it has not been strongly constrained the 
inclination angle for the Crab pulsar and PSR B0540-69. Therefore, we treat the inclina- 
tion angle as a model parameter. The viewing angles ^ of the observer measured from the 
rotational axis is also a model parameter. For this local model in the magnetosphere, the 
current should be dealt as a model parameter and the position of the inner boundary de- 
pends on the assumed current (Hirotani et al. 2003; Takata et al. 2004, 2006; Hirotani 
2006). Instead of the current, however, the ratio of the radial distance to the inner boundary 
and distance to the the null surface r„ is parameterized and is assumed to be constant for 
each field line, that is, rj„(0)/r„(0)=constant, such that the inner boundary locates far away 
from the stellar surface if the null surface locates far away. For example, from Figure [H the 
outer gap accelerator with the non-dimensional current j ~ 0.22 has the inner boundary 
at rin{(j))/r^{(l)) ~ 0.65 in the meridional plane. The altitude of the emission region of the 
secondary pairs is also model parameter, because the magnetic field will be modified by rota- 
tional and plasma effects near the light cylinder and because the last-open field lines may be 
different from the traditional magnetic file lines that are tangent to the light cylinder in the 
vacuum case. To specify the upper surface of the outer gap, it is convenient to refer to the 
footpoint of the magnetic surface on the star and to parameterize the fractional polar angle 
(^f = du/dici where 9u and Oic are the polar angle of the footpoints of the magnetic surfaces 
for the gap upper surface and the last-open field line in the vacuum case, respectively. We 
constrain the boundary of the radial distance to the emission region to r = i?;^- In this 
paper, we apply the model to the Crab pulsar in section Wl\ and to a Crab-like pulsar, PSR 
0540-69, in section i^l 
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4. Model results 

4.1. The Crab pulsar 

For the Crab pulsar, we adopt the inclination angle of a = 50°, the viewing angle of 
^ = 100° and the position of the inner boundary described by rj„(0)/r„((/)) = 0.67, which 
were chosen in Takata et al. (2007). In this paper we chose the fractional angle of a/ = 1 
to explain the phase-resolved spectra, that is, we assume the gap upper surface with the 
magnetic field lines that are the conventional last-open field lines in vacuum. The opening 
angle of the active outer gap in the azimuthal direction is set at 50 = 250°, which as assume 
in section [STTl 

Figure [3] is the photon mapping of the outwardly propagating photons, where the emis- 
sion direction tangent to the local field lines, which is described by = 1, were temporary 
assumed. Figure H] shows the variations of the typical radial distance to the emission points 
of the photons measured by the observer with the viewing angle ^ = 100°. The dashed-line 
in Figure H] represents the radial distance to the emission points that locate beyond the null 
surface on the magnetic field lines coming from the north pole. And, the dotted-lines show 
the distance to the points that locate below the null surface on the magnetic field from the 
south pole. In the traditional study, only the emission beyond the null surface (dashed-fine) 
have been considered. We will see that the emission component below null surface is required 
to explain the phase-resolved spectrum of Peak 1 (Figure [8]). To calculate the phase- resolved 
spectra, we define the phase intervals of Peak 1, Bridge, and Peak 2 as 0.06 — 0.16, 0.29 — 0.4, 
and 0.49 - 0.6 (Figure IH). 



4- 1.1. Pulse profile and Polarization 

Figures [5] and [6] show the predicted variations of the intensity (upper), the position angle 
of the electric vector of the radiation (middle) and the degree of the polarization (lower) as a 
function of the pulse phase from optical to 7-ray bands. To compare with the observed pulse 
profiles in Kuiper et al. (2001), the results were calculated by integrating the photons within 
the energy interval 1-10 eV, 0.1 - 2.4 keV, 20 - 100 keV, 100 - 315 keV, 0.75 - 10 MeV and 
30 — 100 MeV. In the figure, we define the rotation phase $ = in abscissa axis as lying the 
south pole, and the zero degree in the position angle of the bottom panel is corresponding 
to the direction of the rotation axis projected on the sky. 

In Figures [5] and [6], we see that the calculated pulse profile morphology changes signifi- 
cantly as a function of the photon energy likewise the observational pulse profiles (see figure 5 
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in Kuiper et al. 2001). One can see in the results that the Peak 1 emissions dominate in the 
pulse profile in optical bands (left column in Figure [5]). But, the intensity ratio of Peak 1 and 
Peak 2 decreases with increase of the photon energy in optical to soft X-ray bands (middle 
and right columns in Figure E]), and the ratio becomes unity around 100 keV. Eventually, 
the Peak 2 emissions dominate the Peak 1 emissions in hard X-ray and soft 7-ray bands as 
left and middle columns of Figure [6] show. In the hard 7-ray bands, furthermore. Peak 1 
again dominates in the pulse profile (right column in Figure [6]). This predicted morphology 
change of the pulse profiles reproduces the observation very well. 

With the present outer gap model, the morphology change of the pulse profiles is ex- 
plained by the radiation properties beyond the null surface, because the major part of the 
total photons are produced there. The intensity of Peak 1 is much stronger than that of 
Peak 2 in optical wavelengths, because pileup of the photons emitted from the various posi- 
tions due to the special relativistic corrections (i.e. the aberration of the emitting direction 
and the flight-time) occur more efficiently at the Peak 1 phase than the Peak 2 phase. This 
effect has been considered in the previous studies (Romani & Yadigaroglu 1995; CRZOO) 
with a constant emissivity. But, the present calculation produces a Peak 2 emission, which 
becomes more and more pronounced with increasing energies, and eventually dominates the 
Peak 1 emission in the hard X-ray and the soft 7-ray bands. This is because many photons 
of Peak 2 are emitted near the stellar surface and on the other hand, most of the photons for 
Peak 1 are emitted near the light cylinder, as dashed-line in Figure H] shows. As a result, the 
phase-resolved spectrum of Peak 2 of the synchrotron photons is harder than that of Peak 1 
so that Peak 2 dominates in the pulse profile in hard X-ray and soft 7-ray bands. And 
finally, the Peak 1 emissions dominate in the hard 7-ray region, where the inverse- Compton 
process is major emission process (Figure [8]). For the inverse Compton scattering with the 
synchrotron background photons, the averaged strength of the magnetic field affects the 
power of radiation following equation (jl]). Because the averaged magnetic field strength in 
the emission region for Peak 1 is smaller than that for Peak 2, the power of the inverse 
Compton scattering in the emission region in Peak 1 is larger than that in Peak 2, and as 
a result Peak 1 dominates in the pulse profile in 30-100 MeV bands. Thus, we find that 
the synchrotron radiation and inverse Compton scattering of the pairs with the outer gap 
accelerator model can naturally explain the complex morphology change of the observed 
pulse profiles. 

We see a small peak feature at the leading phase of Peak 1 of the pulse profile of 0.75- 
10 MeV bands (middle column in Figure E]). We will discuss this feature in section O together 
with the dependency of the results on the model parameters. In fact, this small peak consists 
of the emissions below the null surface. 
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For the polarization characteristics, Takata et al. (2007) discussed the predicted charac- 
teristics in optical bands (left column in Figure [5]) are in general consistent with the observed 
features (Kanbach et al. 2005) such as a large position angle swing of the electric vector of 
the radiation at each peak and ~ 10% of the degree of the polarization between two peaks. 
In the present study, we find that a similar pattern of the polarization position angle (middle 
panels of Figures [5] and [6]) is predicted in the higher energy bands where the synchrotron 
process dominates. But, we find that the degree of the polarization in the Bridge emissions 
depends on the photon energy. As Figure shows, the polarization degree in the Bridge 
emissions decreases from 10% of the optical wavelengths with increase of the photon en- 
ergy. This is because the tertiary pairs, which have a small Lorentz factor, contribute to the 
bridge emissions with about 10 % of the degree of the polarization in the optical wavelengths 
and contribute with a smaller emissivity in the higher energy region. In 10-100 MeV bands 
(right column in Figure [H]), where the inverse Compton scattering dominates, the degree of 
the polarization is about 10% is predicted. By comparing the calculated position angles 
of the Bridge emissions below IMeV and above 10 MeV, we find that the present model 
predicts that the observed radiations in 10-100 MeV bands are polarized to the direction 
orthogonal to the polarization plane of the radiations below 1 MeV. 

Above 100 MeV, the attenuated curvature photons from the outer gap explain the 
observation as we discussed in section [2l Figure [7] summarizes the polarization characteristic 
of the curvature radiation, where we take into account the acceleration of the particle motion 
due to the curvature of the field lines and of the co-rotating motion for calculating the 
polarization characteristic. We assume 80 % of the polarization degree for each radiation. 
The gyration motion can be ignorable in the outer gap accelerator because the pitch angle 
of the primary pairs becomes very small sin^p ~ 10"''. We can see from Figure [7] that 
the intrinsic level of the degree of the polarization at the radiating point is preserved in 
the Bridge phase for the curvature radiation and a small depolarization appears around the 
peaks. This less depolarization is because the curvature radiation does not related with the 
gyration motion of the pairs. This characteristics of the polarization degree of the curvature 
radiation were also predicted by the caustic model in Dyks et al (2004). The present model 
predicts that the radiation above 100 MeV of the curvature radiation is highly polarized 
more than the radiation below 100 MeV of the synchrotron and inverse-Compton process. 



4-1.2. Phase-resolved spectra 



Figure [8] compares the predicted and observed phase-resolved spectra of the emissions 
of the Peak 1 (left). Bridge (middle) and Peak 2 (right) phases for the Crab pulsar. The 
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dashed-line and the dotted-line are, respectively, the spectra of the synchrotron emissions of 
the secondary pairs beyond and below the null surface. The dashed-dotted line represents 
the spectrum the inverse Compton scattering and includes the emission both beyond and 
below null surface, and the dashed-dotted-dotted hne represents the appearance curvature 
spectrum of the primary particles in the gap . The solid-line represents the spectrum of the 
total emissions. . To compare the results with the data (Kuiper et al. 2001), we normalize 
the calculated spectrum of Peak 1 to data, but we did not normalize individually on the 
pulse phase. 

From the left panel in Figure [H we see that the synchrotron emissions below (dotted- 
line) and beyond (dashed-line) the null charge surface contribute to the total emissions of 
Peak 1 and Peak 2. Although the flux of the emissions below the null surface is a smaller 
than that of the emissions beyond the null surface, we find from the figure that the emissions 
below null surface push up the total flux between 1 MeV and 10 MeV bands for Peak 1 so 
that the calculated phase-resolved spectrum is in general consistent with the data over the 
optical to 7-ray bands. If we consider only the emissions beyond the null surface likewise the 
traditional outer gap model, we obtain calculated total flux of Peak 1 that is significantly 
smaller than the observed flux between 1 MeV and 10 MeV bands. This is because only 
emissions around the light cylinder contribute to Peak 1 as the dashed line in Figure H] shows. 
Because the typical energy of the synchrotron emissions is below 1 MeV with B ~ 10^ Gauss 
of the magnetic field around the light cylinder, its flux decreases exponentially above 1 MeV 
as dashed-line in the left panel of Figure [S] shows. On the other hands, we see that the 
synchrotron emissions below the null charge surface contribute to the emission above 1 MeV. 
This is because the emission region below the null surface for Peak 1 spreads from near the 
stellar surface to near the light cylinder and because the synchrotron emissions with higher 
magnetic field near the stellar surface contribute to the spectrum. 

For Peak 2 (right panel in Figure [H]), we find that the observed emissions below 10 MeV 
are explained by only the synchrotron emissions beyond the null surface, unlike with the 
Peak 1 emissions. For Peak 2 phase, the radiating points from near the stellar surface to 
near the light cylinder contribute to the spectrum as the emissions beyond the null surface 
as Figure H] shows. Therefore, the spectrum of the synchrotron photons emitted beyond 
the null charge surface extends above 1 MeV, as dashed line in right panel of Figure [H] 
shows. The emission region below the null surface locates only near the light cylinder so 
that the emissions (dotted-line in the right column of Figure |8]) is negligible compared with 
the emissions beyond the null surface. 

In the Bridge phase, one can see in Figures H] and [8] that the photons emitted only 
beyond the null surface are measured by the observer with the viewing angle ^ = 100°. This 
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is due to the geometry of the dipole magnetic field, and this feature does not depend on the 
position of the location of the inner boundary, providing that we consider only the outward 
radiations. Namely, the outward emissions contribute to only the Peak 1, Peak 2 and off- 
pulse emissions. Therefore, the present model predicts that the Peak 1 and Peak 2 emissions 
consist of the two components (i.e. emissions below and beyond null surface), while the 
Bridge emissions consist of only one component (i.e. emissions beyond null surface). This 
prediction may be consistent with the observational spectral study done by Massaro et al. 
(2006), who fitted the phase-resolved spectra below 10 MeV with specific two functions for 
Peak 1 and Pea 2, while one function for Bridge. 

When we normalize the model phase- resolved spectra to the observation of Peak 1, we 
find that the calculated fiux of the synchrotron emission in Bridge becomes large compared 
with the observation as middle panel in Figure M shows. We consider that this problem will 
be related with the three-dimensional structure of the emission region (section [5]). 

Between 10 MeV and 100 MeV bands, the emissions from the inverse- Compton scatter- 
ing explain the observation. The present model predicts the fiux ratio of 1 MeV emissions 
and 100 MeV emissions in Bridge phase is larger than those in Peak 1 and Peak 2. This 
feature is consistent with the observations. With the present synchrotron emission and the 
inverse- Compton scattering model, the fiux ratio of 1 MeV emissions and the 100 MeV emis- 
sions is proportional to square of the strength of the magnetic field. The average magnetic 
field in the emission region for Bridge is larger than those for Peak 1 and Peak 2, because 
the emission region of Bridge locates nearer the stellar surface than those of the Peak 1 and 
Peak 2 (Figure Hj). As a result, the predicted fiux ratio is larger in the Bridge phase than 
the Peak 1 and Peak 2 phases. 

From Figure [HI above 500 MeV bands, the attenuated curvature radiation from the outer 
gap (dashed-dotted-dotted-line) explains observation. 

4.2. PSR B0540-69 

PSR B0540-69 is sometimes referred to as a twin of the Crab pulsar. Its rotational 
period, surface magnetic field, and spin down luminosity are similar to those of the Crab 
pulsar. Both pulsar have similar spectral features in X-ray bands. The observed fiux ratio (< 
0.01) of PSR B0540-69 to the Crab pulsar will represent ratio of the distances to the pulsars, 
which are c? ~ 49 kpc for PSR B0540-69 and ci ~ 2 kpc for the Crab pulsar. Furthermore, 
the radio giant emissions, which may be associate with the high-energy emissions, have been 
observed from both pulsars (Johnston et al. 2004). While they have theses similarities, the 
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observed pulse profiles are quite different to each other. The pulse profile of PSR B0540-69 
consists of a single broad peak and top of the pulse shape becomes relatively flat, whereas 
the pulse profile of the Crab pulsar shows a sharp double-peak structure, de Plaa et al (2003) 
suggested that the observed pulse profile of PSR B0540-69 is result of two pulses separated 
0.2 in phase. These difference will be useful to examine the model capability. 

The magnetospheric X-ray photon number density will be proportional to rix oc Lgap/Rf^, 
where Lgap ~ f{Ric/2)^Lsd and Lg^ ~ 3.85 x lO^^P'^B"^ erg/s is the spin down energy. The 
ratio of the value of the mean free path, A oc l/ux, of a primary curvature photon in the 
magnetospheres of PSR B0540-69 to the Crab pulsar is estimated as ~ 2. Because the 
thickness of the emission region of the secondary pairs is expected to be related with the 
mean free path of the pair-creation, PSR B0540-69 will have a thicker emission regions of 
the secondary pairs than the Crab pulsar. 

Zhang & Cheng (2000) set the gap thickness of PSR B0540-69 as A ~ 0.21Ric, and they 
chose the inclination angle a ~ 50° and the viewing angle ^ = 76°. The dashed-line in left 
panel of Figure [3 shows the pulse profile in 10-20 keV bands for the inclination angle a ~ 50° 
and for the viewing angle angle ^ = 76° and ^ = 104°. The present model produces a broad 
single peak in the pulse profile as Zhang & Cheng (2000) showed. In this viewing angle, 
however, we see that the fiat feature of the top of the pulse profile appear with a shorter 
width (~ 0.1 in rotational phase) than the observational value ~ 0.2. A viewing angle more 
far from 90° produces a more narrow pulse. On the contrary, if we adopt a closer viewing 
angle to 90°, the produced pulse profile has two distinct sharp peaks. 

With the dipole magnetic field configuration, a more consistent results with the ob- 
servation is obtained with a smaller inclination angle than a = 50°. The solid-line of the 
left panel in Figure M shows the pulse profile of 10-20 keV bands with the inclination angle 
a = 30° and the viewing angle ^ = 98.5° (and ^ = 81.5°). We see that a wider pulse width 
is produced and the shape of the top of the pulse profile is more similar with the observed 
profile (see figure 2 of de Plaa et al. 2003). Therefore, the present model predicts that 
PSR B0540-69 has a thicker emission region (about factor two) than the Crab pulsar and a 
smaller inclination angle than a = 50°. 

The expected phase-averaged spectrum in optical to 7-ray bands is displayed in right 
panel of Figure [9] and is consistent with the observations except for the UV data. Figure [TO] 
shows the calculated pulse profile in the optical bands. We see that the pulse width of 
the optical bands becomes wider than that of the X-ray bands and the narrower first peak 
becomes more stronger than the narrower second peak. These features are also consistent 
with the data. We also display the predicted polarization characteristics in the optical bands 
for ^ = 98.5° in Figure [TOl For the viewing angle ^ = 81.5°, which is mutually symmetric 
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with respect to the rotational equator with ^ = 98.5°, the expected pulse profile and the 
degree of the polarization are identical with those in Figure [10], but the position angle of the 
polarization is mirror symmetry of that in Figure [TU] with respect to the equator 0°. 

5. SUMMARY AND DISCUSSION 

In summary, we studied non-thermal emission process of the Crab pulsar with the 
two-dimensional electrodynamical and the three-dimensional studies. We calculated the 
synchrotron and the curvature radiation process of the primary particles in the gap using 
the electrodynamical study and showed a result which produces a consistent gamma-ray 
spectrum with EGRET observation. Base on the two-dimensional electrodynamical study, 
we conducted the three-dimensional model to calculate the synchrotron and the inverse 
Compton scattering of the secondary pairs. We computed the pulse profiles, spectra and 
polarization characteristics in optical to 7-ray bands to compare the data of the Crab pulsar 
and PSR B0540-69. In the three-dimensional study, we demonstrated that the synchrotron 
radiation and the inverse Compton scattering with the outer gap geometry naturally explain 
the morphology change of the observed pulse profiles for the Crab pulsar. We also showed 
that the observed phase-resolved spectra of the Crab pulsar are explained by the synchrotron 
radiation and the inverse Compton scattering of the pairs below and beyond the null surface. 
The predicted polarization characteristics of the Crab pulsar were in general consistent with 
the data in the optical wavelengths. We also found that a broad single pulse profiles and 
phase-averaged spectrum of PSR B0540-69 are explained by the present outer gap model 
with a thicker emission region and a smaller inclination angle than the Crab pulsar. 

We have shown that the present outer gap model can explain a lot of the observational 
properties. But, a further discussion for the emission of soft 7-ray regions will be required, 
because we found that the calculated pulse profiles of soft 7-ray regions, say 0.75-10 MeV, 
easily has a triple-peak structure in a single period, as the calculated pulse profile of 0.75- 
10 MeV in Figure E] indicates. In this section, therefore, we discuss the pulse profile of the 
0.75-10 MeV bands together with the dependency of the results on the model parameters. 

As Figure H] shows, the both emissions beyond and below the null surface contribute to 
Peak 1. In this case, the small leading peak, which leads the mail first peak, easily appears 
in the pulse profile of 0.75-10 MeV, because the synchrotron spectrum below the null surface 
is harder than that below the null surface. Below MeV bands, the intensity of the small peak 
is much weaker than that of the main first peak, and therefore it can not be identified in the 
pulse profile. However, the emissions beyond null surface, which are major contribution of 
the Peak 1 emissions, are soft so that the intensity decreases exponentially above 1 MeV as 
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the calculated phase- resolved spectrum (the dashed line in left panel of the Figure [8]) shows. 
On the other hands, the spectrum of the emissions below the null surface extends above 
1 MeV because the photons were emitted near the stellar surface. As a result, the peak by 
the synchrotron emissions below the null surface becomes conspicuous in the pulse profile 
above 1 MeV. 

With a viewing angle more close to 90° than the present ^ = 100°, one can obtain a pulse 
profile with only two peaks in over the energy bands. In this case, the leading small peak 
appears the same phase with and are not discriminated from Peak 1 even in the pulse profile 
of 0.75-10 MeV bands. In such a viewing angle, however, the phase separation between the 
main two peaks (Peak 1 and Peak 2) is equal to or larger than 0.5 phase. For a viewing 
angle more far from 90°, on the other hands, the triple-peak structure appears not only in 
the pulse profile of 0.75-10 MeV bands, but also the pulse profiles in lower energy bands. 
On these ground, we consider the viewing angle ^ ~ 100° is better fit parameter for the 
inclination angle of a = 50°. 

The altitude of the magnetic surface for the lower boundary of the pair emission region 
(i.e. the upper boundary of the outer gap) also affects to the pulse profile. When we adopt a 
more higher altitude as the position of the boundary than the present position (i.e. when we 
adopt a smaller value oi aj = QujQic than the present a/ = 1), the phase separation between 
the leading small peak and the main first peak becomes wider. As a result, the pulse profiles 
have three peaks in a single period from optical to 7-ray bands. On the contrary, when we 
adopt a lower altitude for the position of the boundary with a fraction angle aj larger than 
unity, we obtain the pulse profiles having only two main peaks, but the phase separation 
between two-peaks becomes to be equal to or larger than 0.5 phase. As the lower boundary 
of the secondary emission region, therefore, we conclude that a magnetic surface close to the 
surface of the conventional last-opened field lines, a ~ 1, is preferred for the lower boundary 
of the emission region. We also examined the different inclination angle. For example, 
the inclination angle ^ = 60°, the viewing angle C, ~ 100° and aj = 0.99 reproduced the 
consistent phase-resolved spectra and the pulse profiles with the Crab data. However, we 
also obtain a small leading peak in the pulse profile of the 0.75-10 MeV. 

Because the phase of the peak in the pulse profile depends on the configuration of the 
magnetic field, the present results would indicate that the actual magnetic field in the pulsar 
magnetosphere is modified on important level from the vacuum dipole by plasma effects 
(Muslimov & Harding 2005). To examined this issue, we need to solve the magnetic field 
configuration with the current. 

A more quantitative arguments in the three-dimensional structure of the magnetosphere 
will be required to resolve also issue that the calculated Bridge emissions are large relatively 
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to the observed emissions in the soft 7-ray bands, when we normahze the calculated flux 
using the Peak 1 emissions, as Figure [8] shows. This will indicate that the dependence of the 
trans-field thickness of the emission region is described as a more complex function of the 
azimuth than the present model, in which the emission region was assumed to be located 
between two magnetic surfaces. In the actual magnetosphere, the trans-field thickness around 
the meridional plane, where the Bridge emissions take place, may be small relative to the 
thickness of the leading and the trailing sides, which produce the Peak 1 and Peak 2 emissions. 
This will be expected because the null charge surface is located nearer the stellar surface 
around the meridional plane compared with the leading and the trailing side. Therefore, the 
surface X-rays are more dense around the null surface in the meridional plane so that the 
mean-free path of the pair-creation and the resultant trans-field thickness of the outer gap 
and the emission regions are shorter around the meridional plane than that of the leading 
and trailing sides. Since the three-dimensional structures of the acceleration region and the 
emission region have not been solved up to now, the above issues will demand us to perform 
a more quantitatively study of the three-dimensional structure in the subsequent papers. 

Finally, we note that it is important to observationally constrain the inclination angle. 
Takata et al. (2007) searched the parameter range to explain the polarization characteristics 
and the pulse profile in the optical wavelengths of the Crab pulsar. In their study, a wide 
range of the viewing angle for each inclination angle was allowed for the possible parameters 
to explain the optical pulse profile and the polarization. In this study, we found that the 
narrower range of the viewing angle can produce the consistent pulse profile and phase- 
resolved spectra with data in optical to 7-ray bands. Therefore, if the inclination angle can 
be determined in other ways, the viewing angle will be strongly constrained by the present 
study. 

Authors appreciate fruitful discussion with K.S. Cheng, K.Hirotani, S.Shibata and 
R.Taam. Authors also thank anomarous referee for insightful comments on the manuscript. 
This work was supported by the Theoretical Institute for Advanced Research in Astrophysics 
(TIARA) operated under Academia Sinica and the National Science Council Excellence 
Projects program in Taiwan administered through grant number NSC 96-2752-M-007-001- 
PAE. 



A. Stokes parameter of the inverse Compton scattering 

In this appendix, we derive the Stokes parameters given by equations ([8]) and ffTOj) . which 
are respectively the Stokes parameters of the inverse Compton process for the unpolarized 
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and the polarized background radiation. 

In this appendix, as polar-axis (z-axis), we choose the direction of the particle motion 
at local point. With equation ([3]), the unit vector e* along the z-axis is described by the 
unit vectors (without asterisk) in the coordinate system based on the rotation axis as 

e* = po cos 6pb + po sin ^p6_L + /3coe<^, (Al) 

where 6p is the pitch angle, b is the unit vector along the magnetic field, and b± = 
±(cos 6(j)K+ sin 6(j)K x 6) is the unit vector perpendicular to the magnetic field line, K is the 
unit vector of the curvature of the magnetic field, and 6(f) is the phase of the gyration motion. 
The parameter /3o is determined from |e*| = 1. W define x and y axes with (Figure [TT!) 

el = [(3oK - {l3,Kjhr)er]/\[l3,K - (/5,ir<^/6,)e,]|, and = e: x el (A2) 

The propagating direction of the background radiation, fco, are represented by a polar angle 
6q and a azimuthal angle 0O) that is, = sin 6^0 cos 0oe* + sin^o sin(/)oe* + cos 6*06*. The 
propagating direction of the scattered radiation, ki, is represented by coordinates (6*1, 0i). 
If a particle moves with the Lorentz factor Fe, we relate the quantities of the rest frame of 
the electron (with prime) and of the observer frame (without prim) as 

cos ^' = (cos ^ -/?)/(!-/? cos ^), 0' = 0, 

for the polar and azimuthal angles and e = F(l + (3cosd')e' for the photon energy. If 
the background radiation is polarized, we describe the direction of the polarization with 
Eyj = sin 9u, cos 0«,e* + sin 9^, sin 0^e* + cos 6'^^,e*. 

In the calculation, we first compute the Stokes parameters of the Compton process in 
the electron rest frame, and then we perform the Lorentz transformation to the observer 
frame. 



A.l. Unpolarized background field 

For the Compton scattering with the unpolarized background radiation, the Stokes 
parameters, which are represented in the coordinate basis orthogonal to and parallel to the 
scattering plane, in the electron rest frame are 



Pi 

u' 

P\\ 



-T + -r — sm w\ 
sin^ w' 



(A3) 
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(McMaster, 1961), where Sq corresponds to the emissivity, p^' and describe the po- 
larization direction orthogonal to and the parallel to the plane of the scattering, respec- 
tively. The dimensionless quantities e'^ and e[ are, respectively, the energy of the background 
and the scattered photons in units of the electron rest mass energy and are connected by 
e[ = ey[l + €0(1 — cosWg)]. The scattering angle w'g is defined by 

cos w'^ = feg ■ k[ = sin 9'q sin 9[ cos(0o ~ 4>i) + cos 9^ cos 9[. (A4) 

For the unpolarized background radiation, we find that the scattered radiation is partially 
polarized orthogonal to the scattering plane. 

We transform the back ground frame to the one related with the particle motion direction 
(z-axis) projected onto the sky. Using such coordinate system, the Stokes parameters become 

Qu = p1 cos 2ri' — sin 21]' = sin^ w'^ cos 2//', (A5) 

and 

u'^ = p\ sin 27]' + p\ cos 2'q' = sin^ w'^ sin 2//', (A6) 

where t]' is angle between the projected z-axis on the sky and the direction orthogonal to the 
scattering plane, and its cosine is given by cos?]' = — sin 6*3 sin(0Q — (f)[)/ sinw'^ (Figure [TTi). 

Writing the unpolarized background radiation in the rest frame with /^'(ecfep), the 
Stokes parameters of the total amount of the scattered radiation are give by i^'(e[,k[) = 
(3aT/167r)//,"'s^'rffi', Q'^{e[,k[) = {3aT/lQ7r) J IJ}' q'jn' , andK{e',,k[) = {3aT / IQtt) J I^' u'JQ' , 
respectively. If we take the Thomson limit, eg = e[, we arrive the equations (2), (4a) and 
(4b) of Begelman and Sikora (1987). 

By performing the Lorentz transformation from the electron rest frame to observer 
inertial frame, and by transforming the background coordinate to one related the projected 
rotation axis on the sky, the Stokes parameters become 

^ik^,e^) = Dlt''{k[,e,) 

q{k,,e,) = Dl[q'{k[,e[)cos2C ~u'{k[,e[)sm2C], (A7) 
«(fei,ei) = Dl[u'{k[,e[)cos2C + q'ik[,e[)sm2C], 

where D2 = ei/e'i = r^^(l — /5cos6'o)~^ is the Doppler factor, and the azimuth angle ( is 
defined by angle between the angle between the direction of the rotation axis and of the 
particle motion projected on the sky. 

Finally, integrating i, q and u over the particle momentum distribution, the Stokes 
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parameters are given by 



dlu{ki,ei)/ds 
dQu{ki, ei)/ds 
dUu{ki, ei)/ds 



dVp 



drie 
dT, 



i(fci,ei), 

;i -/?cos^o) X { q{ki,ei), 
u{ki,ei), 



(A8) 



where the extra factor (1 — /? cos ^o) represents relative motion of photons and particles along 
the photon motion. 

Now, we suppose the total background radiation in the observer frame that Ib{ko, eo) = 
CoeQ"6{9 — 9o)6{(j)—(j)o)- In the rest frame of an electron, the background radiation distributes 

as 

i3w> ^_ C'oe[r" 



Il{k'„e',) = Dihiko,eo) 



(A9) 



ri+"(i + ^cos^')^" 

where we use = e'^/eo = T-\l+(3 cose')-\ and 6 {6 -6 0)6 {(f)- (po) = 

which is because 6(6 - eo)6{(f) - 0o)rffi = 5(6' - ^o)<^(0' - (l>o)d^' and dil = Djdn'. 

By denoting the polarization degree of the total background radiation as Hgyn, that is, 
if = {1 — Ilsyn)Ib^ finally we arrive equitation (IHl), 



dlu{ki, ei)/ds 
dQu{ki,ei)/ds 
dUu{ki, ei)/ds 



n)^Co / dV, 
i67r 



dUe 
dVf. 



9Ufc'i>e'i)cos2C - <(fci,ei)sin2C, 
g;(fc;,e;)sin2C + <(fci,ei)cos2C. 



r4+°(l -/? cos + cos 0^)1+" 

(AlO) 



A. 2. Polarized background field 



For the polarization background radiation, the Stokes parameters corresponding to equa- 
tion (1A3I) become 

4i + 4 - sin^ e[ cos^ a; 



sin' e[ - (1 + cos2 e;)cos2A; 



P\\ ) 2cos^isin2A;, 

where A^ is the angle between the polarization plane of the background radiation 
orthogonal plane to the scattering plane and its cosine is given by 



(All) 
and the 



cosAp = (sinu;^) ""^ [sin 6'q cos 6*^ sin sin(0o — </>«)) + cos sin ^'^ sin 6^^ sin((/)^ — 
+ sin 9'q sin 9[ cos 9'^ sin(0i — 0o) , 



(A12) 
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where 9'^ and 0^ represent the polar angle and the azimuthal angle of the polarization plane 
of the background radiation, respectively. 

Transforming the background coordinate to one related with the particle motion, the 
components of the Stokes parameters representing the polarization are, respectively, given 
by 

= [sin^ w[ — {1 + cos^ w') cos 2Ap] cos 2-r]' — 2 cos' w sin 2Ap sin 2r], (A13) 

and 

u'p = [sin^ w[ — {1 + cos^ w') cos 2Ap] sin 2t]' + 2 cos' w sin 2Ap cos 2r]. (A14) 

After same procedure with the case of the unpolarized background radiation (or replacing 
subscript u and 1 — Hsyn in the equation flAlOp to p and to Hsyn), we arrive equation ffTOj) . 
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Radial distance [Light radius] 



Fig. 1. — Electric structure in the gap. The sohd hne shows the solved accelerating field 
along the field hne, which locate at 50 % of the thickness from the last open field hne. This 
results are for The inclination angle is 50° and for about 22 % of the Goldreich- Juhan current 
in the gap. Beyond the null charge point, the dashed line is vacuum electric field obtained 
by Cheng et al (1986a). 
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Fig. 2. — The spectrum of the synchrotron and the curvature processes of the primary 
particles in the outer gap. The sohd hne represent the intrinsic flux from the outer gap. 
The dashed line show the appearance flux after pair-creation process outside of the outer 
gap with the soft photons emitted by the synchrotron radiation of the secondary pairs. The 
observation data are taken from Kuiper et al. (2002) and reference therein. 
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Fig. 3. — Emission region projected onto the $)-plane for the magnetic surface af = 1. 
Here,the emission direction is tangent to the direction of the local magnetic field. The 
inclination angle is a = 50°, the emission region extends from rj„ = 0.67 to r = Ric, and the 
width of the polar cap angle is about 250°. The circles show the shape of the polar cap. The 
viewing angle ^ = 100° (horizontal dashed- line) is used in the calculations for Figure Hh^lHl 
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Fig. 4. — Variation of radial distance with pulse phase for the inchnation angle a — 50° 
and the viewing angle ^ = 100°. The dashed-hne and the dotted-line represent the emission 
regions beyond and below the null surface, respectively. 
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Fig. 5. — Variation of the intensity (tot), the position angle (middle dots) and the degree of 
the polarization (bottom solid-line) as a function of the rotational phase. For the reference, 
the light curve is overplotted at middle and bottom panels (thin dashed-lines). To compare 
with the observe pulse profiles of Kuiper et al. (2001), the results were obtained by integrat- 
ing the photons within the energy interval 1 — 10 eV (left column), 0.1 — 2.4 keV (middle 
column) and 20 — 100 keV (right column). The model parameters are a = 50°, ^ = 100°, 
rin — 0.67r„ and af — 1. 




Fig. 6. — Same with the Figure O but the resuhs are for the interval of photons energy, 
100-315 keV (left), 0.75-10 MeV (middle) and 30-100 MeV (right) 




Fig. 7. — The polarization characteristic of the curvature radiation in the outer gap. The 
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Fig. 8. — Phase-resolved spectra for Peak 1 (left), Bridge (middle) and Peak 2 (right) for 
the Crab pulsar. The dashed-line and the dotted-Iine show the spectra of the synchrotron 
emissions of the pairs beyond and below the null surface, respectively, the dashed-dotted-line 
show the spectrum of the inverse Compton scattering, and the dashed-dotted-dotted-hne 
show the attenuated curvature spectrum in the outer gap. The solid-line represents the 
spectrum of the total emissions. The model parameters are a. — 50°, ^ = 100°, rj„ = 0.67rn 
and aj — \. The observation data are taken from Kuiper et al. (2001). 
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Fig. 9. — Left; The pulse profile in 10-20 keV bands for PSR B0540-69. The solid-line 
is result for a = 30° and ^ = 98.5° (and ^ = 81.5°), and the dashed-line is for a = 50° 
and ^ = 104° (and ^ = 76°). Right; The phase-averaged spectrum of PSR B0540-69 for 
a — 30° and ^ = 98.5° (and ^ = 81.5°). The solid-line represents the total emissions of the 
synchrotron radiation (dotted-hne) and the inverse-Compton scattering (dotted-line) . The 
observation data are take from de Plaa et al (2003). 
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Fig. 10. — Pulse profile in optical bads of PSR B0540-69. The result is for a = 30° and 
^ = 98.5°. The top, middle and bottom panels show the intensity, the position angle of the 
polarization and the degree of the polarization, respectively. 
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Fig. 11. — Left: Coordinate system used to describe scattering in the electron rest frame. 
The polar axis, e*, is chosen to coincide with direction of the electron motion in the observer 
frame. The definition of e* and e* are seen in equation flA2p . The vectors feo and ki are the 
directions of the background and the scattered radiations, respectively. If the background 
radiation is polarized, the vector E^] is used to represent the polarization plane of it. Right: 
Coordinate system on the sky. VLp and e*^ is, respectively, the directions of the rotational 
axis and of the particle motion projected on the sky. The position angle of the polarization 
plane (Egw) is measured anticlockwise from the direction of the rotational axis on the sky. 



